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INFLUENCE OF MAGNETIC SHEAR ON THE COLLISIONAL CURRENT DRIVEN 
ION CYCLOTRON INSTABILITY 


I. INTRODUCTION 

Intense field aligned currents observed In the auroral E and F regions 
(Park and Cloutier, 1971, Cahill et al., 1974, Kelley et al., 1975) were 
conjectured to be responsible for the generation of electrostatic Ion 
cyclotron waves. Since the Ionosphere Is colllsional in the E and F regions 
the field aligned currents give rise to negative energy waves which grow due 
to the dissipation associated with electron-neutral collisions (Chaturvedl, 
1976). These are electrostatic waves and lead to density irregularities in 
the medium. In fact, radar measurements in the E region detect 
irregularities collocated with auroral electrojets (Greenwald et al., 1975), 
and a.c. field measurements (Kelley et al., 1975) have been identified as 
being due to unstable ion cyclotron waves (Drummond and Rosenbluth, 1962 - 
collisionless domain; Kindel and Kennel, 1971 - weakly colllsional 

domain). Sounding rockets launched from Syowa station, detected 
transversely propagating electrostatic plasma instabilities of scale sizes 
around 200 km in association with strong field aligned currents (Ogawa et 
al., 1981). Chaturvedl (1976), using a local fluid analysis, showed that, 
in a partially ionized plasma, field aligned currents can support almost 
transversely propagating ion cyclotron waves owing to resistivity 
experienced by electron parallel motion (electron-neutral collisions), while 
the ion-neutral collisions were found to have a mild stabilizing effect. 
The field aligned currents observed are usually ** 1-10 pA/m (Kelley et al., 

1975), although recently Burke et al. (1983) reported observations of large 
2 

currents - 135 pA/ra . These field aligned currents can generate magnetic 
shear which in turn affect the mode structure. This is especially so in the 
case of collisionless current driven ion cyclotron waves (Ganguli and 
Bakshi, 1982). 

In this brief report we investigate the effects of a self-consistent 
magnetic shear on the colllsional current driven ion cyclotron instability 
(CIC). We find that the magnetic shear has no noticable stabilizing effects 
on the current driven ion cyclotron instability in the domains of interest 

pertaining to auroral conditions (characterized by V, - 0.5 - 5 km/s, 

2 Q 4-1 

corresponding to currents of order 10 - 100 pA/m , and v ^ 10 s , where 
Manuscript approved March 20, 1984. 
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is the electron drift velocity along the magnetic field and is the 

electron-neutral collision frequency). This is because the self-consistent 

magnetic shear produced by the currents observed in the auroral region is 

very small. This result is in contrast to the major effects found by 

Ganguli and Bakshi (1982) in the collisionless domain. Their kinetic 

treatment shows that, due to the nonlocal boundary conditions, magnetic 

shear affects the collisionless current driven ion cyclotron instability in 

two ways: One, due to an explicit shear dependent term of 0(p /L ) which 

s 3 

vanishes for P s /L g + 0* Two, due to a shear independent term containing the 

derivatives of the potential function defined in equation 18. In the fluid 

example, the second term happens to be zero (since the derivative of Q, 

defined in Eq. 18, vanishes at Xq which is defined in Eq. 19) and hence the 

lack of the significant effects in the p /L ♦ 0 limit. Furthermore, 

s s 

for p /L « 1 magnetic shear localizes the wave packet around k /k « 1, 
s s z y 

Indicating that the mode is almost perpendicularly propagating, and stronger 
magnetic shear stabilizes the instability. The paper is organized as 
follows. In the next section we give the theory and derive the nonlocal 
mode structure equation. In the third section we present the results, and in 
the last section we apply the results to the auroral ionosphere and discuss 
future work. 


II. THEORY 

We consider a partially ionized plasma such as the one encountered in 
the auroral E region. The geometry used in this analysis is as follows. 
The magnetic field is aligned with the z-direction along which flow the 
currents that drive the CIC. Because of the self-consistent shear produced 
by the field aligned currents, the magnetic field acquires a component along 
the y-direction as x is varied, 

B(x) * B (£ + x/L g . (1) 

This implies that the parallel wave vector becomes a function of x. 


k 

z 


k z + V X/L s J ’ 


( 2 ) 
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where L g is the shear length. The self-consistent shear is calculated from 
the Maxwell's equation 


7x 1-7 L £-“ L t ne V 


(3) 


where is the drift velocity parallel to the z-direction, and is given by 


-1 


L » ■ 0 ' , , | . C S / V 


(4) 


1/2 

where c g * (T^/m^) , p g * c g /fl^, is the ion cyclotron frequency, and 


a = U g /cj 2 


(5) 


where u> g and are the electron plasma and cyclotron frequencies, 
respectively, c is the speed of light, and M and m are the ion (N0 + ) and 
electron masses, respectively. 

The basic equations describing the problem are as follows. The 
electron and ion continuity equations are given by 


3n 


- 0 


3t v ~a-a' 

where a indicates the species. The equation of motion for the ions is. 


( 6 ) 


. V* B 

Mn V. - - Vp. + en( -V$ + —-J - Mnv V, , 

dt —1 i v c ' i —1 


(7) 


while the equation of motion for electrons is. 


V x B 

0 =« - Vp^ - en(-V$ + — —-J - mnv^ , 


( 8 ) 


where v and v are electron-neutral and ion-neutral collision frequencies, 
e l 

respectively, p a » n a ^ a * s c ^ e pressure of species a, and d/dt * 

3/3t + V*7. In this paper we ignore the electron inertia, and assume that 

« l, and « 1. We al30 Ignore the viscosity and thermal 

conductivity terms because these terms only introduce a numerical constant 
that multiplies the growth rate and do not introduce any additional k 2 /ky 
terms (see, for example, Chaturvedi and Kaw, 1975). We also ignore the 
elctromagnetic effects and consider only to electrostatic perturbations. We 
derive the mode structure equation assuming that the perturbed quantities 

3 






vary as f ' f(x) exp [-i(wt - k y - k zJJ . We note that we do not Fourier 

y z 

analyze in the x direction since the magnetic field is assumed to vary in 
that direction. 

From Eq. (8), we obtain the perpendicular and parallel (to B) 
components of the electron velocity 

\ U+v2/fl2j ® 7i * j X ' 

e ej 7 p v < 9 > 

- f ,_ i _£ . £ , _e. 

k mn m Q J 

e 

V - - — [—— - - 7 <j>| . (10) 

ez v L mn m z* J v ' 

e 


From Eq. (10) we define 


V = V° - — (— 7 4 >°\ 
d ez v Vb v z * )> 
e 


From Eq. (7) we obtain for the ions 

i<»>/fl , V 


where we have assumed T g ^ and ♦ — <fr° + <p (<p being the perturbed 
potential). Substituting Eqs. (9) thru (13) in Eq. (6) we obtain, after 
subtracting the electron continuity equation from the ion continuity 
equation, the nonlocal mode structure equation 

. 2 , 2 _2 .2 2 , . 2 . .2 


a 2? k (u -O.-k c j + i(m/Mjk v (ui-k V.J + iuv k^ 
3 ♦ + 1 y s' v ' y e v z d ; _ i z 


k l c t ~ iCm/MJv (a)-k V.J 


where k^ ■ k z + k^x/L^ and where quasi-neutrality has been assumed. The 
local limit of Eq. (14) (L g + » and 3/3x « k y ) yields Q(x) - 0, where Q(x) 


l 


| 






are Che second two terms in the above equation. This gives the results 
obtained previously by Chaturvedi (1976). The linear local growth rate is 

k • V v 

os) 

2M k r 


where 


2 2 2 2 

ur = flj + k c . 
r i y s 


U) = U1 + iY 

r 


Gq. (15) leads to the threshold condition 




(16) 


y e 

Note that the nonlocal treatment allows for the formation of a wave packet 
in the x-direction compared with the plane wave [exp(ik^x)] solutions of 
Chaturvedi (1976). 


We solve Eq. (14) analytically and numerically to determine the 
eigenvalues and eigenfunctions. We follow the techniques developed by 
Ganguli and Bakshi (1982) to solve Eq. (14) as 

^4+Q(x)J-0 f (17) 

3x 



2 

We expand Q(x) around £ = x-x ■ 0 to 0(5 ) where x^ is the value of x for 

o ° 

which the local growth rate is a maximum. 


x 

o 


2 oj 


k 

y d 


(19) 



. • 


i 
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Now, we have 


2 * Q" 

+ U + n 2 J J(n) - 0 (20) 

3n 


where n * 5 + Q'/Q" and a - Q - Q'/2Q". The primes indicate the 
o o o o o r 

derivatives with respect to x evaluated at x ■ x . 

Equation (20), which is in the form of a Weber equation, has the 
solution (for the l * 0 mode) 

<Kx) - $ exp[- y (x - x mo J 2 / x ^i ex P Liq(x)J (21) 


where l is the order of the mode and q(x) is a real function. The 
wavepacket peaks at (Ganguli and Bakshi, 1982) 

x * x - 
mo o 


Re (-Q~/2) 


1/2 


( 22 ) 


The dispersion relation is given by 


Q~ 1/7 (Q ') 2 

Q(x ) - (2 1 + 1) -~J 1/2 P_/L_ + ° 


s' s 2Q"" 
o 


(23) 


It can be easily shown from Eqs. (18) and (19) that Qq( x q ) “ 0. Thus, Eq. 
(23) reduces, for the l « 0 mode, to 


QUJ-l-Q" /2j 1/2 P s /L s . 


(24) 


Equation (23) is anologous to the dispersion relation obtained by Ganguli 

and Bakshi (1982). In Eq. (23) the nonlocal effects are contained in the two 

terms on the right hand side, one of which contains an explicit shear 

dependent term and the other contains a shear independent nonlocal term. 

This equation shows that in the p /L 0 limit the term 

s s 

containing Q' could significantly modify the growth rate. In fact, this 
shear independent nonlocal term drastically reduces the growth rate in the 
collisionless kinetic case. Whereas, Eq. (24) shows that in the fluid case 
any stabilizing infuelnce by the magnetic shear is due only from the shear 
dependent term whiche indicates that large magnetic shear is probably needed 




to stabilize the instability. In fact, the results given in the next 
section support this observation. Equation (24) is solved to obtain the 
eigenfrequencles and the results are presented in the next section. 


III. RESULTS 


We solve the dispersion equation given by Eq. (24) numerically in 
various parameter domains. These results are verified by directly solving 
Eq. (14) using a numerical shooting code. We present the former results in 
the following. In Fig. 1 we plot the growth rate normalized to the ion 
cyclotron frequency (y/0^) versus the normalized wavenumber (k^P g ) for 


Vs 


50, V /SI 
e e 


0.01, and 


0.01. In plotting these results we 
The solid line represents t' 


treated P s /L g as an external parameter, 
shear-free case and the dashed line represents the case where p g /L g - 0.01 
We see from the figure that modes with small k^p g are not strongly affec 
by the magnetic shear. The growth rate for k^p g “ 0*3 reduces from 0.24 ■<* 
0.20 for P g /L g ■ 0.01, whereas the growth rate for k^p g * 0.7 reduces by 
about 30% from 0.75 to 0.55. 


In order to throw some light on the magnetic shear required to 

significantly reduce the growth rate, we plot in Fig. 2, the normalized 

growth rate versus the normalized shear length (p /L ), treating p /L as 

s s s s 

an external parameter. We choose k p *0.3, and 0.5 and plot the growth 

y s 

rate for the parameters of Fig. 1. This figure shows that moderate to 

strong magnetic shears, p /L } 0.01, are required to reduce the growth 

s s 

rate significantly. For example, for k p g - 0.3 the growth rate is reduced 

by about 50% for p /L - 0.04 and the mode is stabilized for p /L > 0.065, 
s s s s 

and for k p » 0.5 the growth rate is reduced by 50% for p /L = 0.04 but 
y s s s 

the mode is stabilized for p /L - 0.078. However, the theory begins to 

s s 

break down for large shear value p /L >0.1. We find that for large 

s s 

shears, the growth rate drops linearly as a function of inverse shear 

-4 

length. For very small shears, p g /L g * 10 the growth rate remains 
essentially constant and the growth rate decreases noticeably for P g /L g > 
0.01. In an ionospheric environment B - 0.5 Gauss and p g - 2.5 meters, and 
Eqs. (1) and (3) show that p g /L g ■ 0.01 and 1.0 correspond to currents of 
0.1 A/m^ and 1 A/m^, respectively. 
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In Fig. 3 we present the wave eigenfunctions for the following 

-2 -2 

parameters: V. ■ 50 c , v /Q ■ 10 , v./Q. - 10 and for k p ■ 0.3, 

d s ei e ’ i i y s 

Y/8. ■ 0.24. The solid line represents the real part of the wave function 

($) and the dashed line |$j. Two important features are to be noted in 

this figure. One feature is that the wavepacket localizes around x Q /p g « 1, 

i.e., for small shear lengths k /k is « 1, indicating that the mode is 

z y 

almost perpendicularly propagating. This can also be seen from Eq. (19) 
which yields 

k /k - 2(io/k V,J(p /L J « 1 for p /L ' 10~ 7 . 
z y y ds s' s s 

The second feature is that the width of the wave packet is of the order 

200 p suggesting a localization region of 500 meters for p * 2.5 meters 

8 + S 
(corresponding to NO ions in the ionosphere). 

The effect of self-consistent magnetic shear is understood by plotting 

the normalized growth rate versus the normalized drift velocity. For this 

purpose, we use Eq. (4) to express the shear length in terms of the drift 

velocity and solve Eq. (24). In figure 4, we plot the growth rate 

for o ■* 10 6 , 10 5 , and 10 * (uj 2 /u> 2 » 10 2 , 10 3 , and 10 4 , respectively; 

pe ce 

c =• 500 m/s) and for v /ft * 0.01, v /ft * 0.01, and k p > 0.3. We also 
3 e e I l y s 

give the shear-free local growth rate (curve A) to compare with the growth 
rate with shear. Three points are to be noted in this figure: (1) a » 10 ^ 
corresponds to ionospheric parameters. We find that the growth rate is not 
too different from the shear-free case. In fact, the growth rate curve 
overlaps curve A and eventually turns around for V^/cg > 500. (2) Curve B, 
which gives the growth rate for a ^ 10 3 , shows that the growth rate drops 

significantly beyond drift velocities > 150 c ; and (3) curve C, which 

-4 

represents the growth rate for a ^ 10 , shows that the optimum drift 

velocity (the drift velocity for which the growth rate is a maximum), 80 c g , 
is much smaller for this case. This leads to the conclusion that the 

optimum drift velocity decreases as the parameter a is increased. At large 

drift velocities, the self-consistent magnetic shear produced is large 

enough that the mode is stabilized and the growth rate tends to zero. 




IV. *■ SCUSSION 


We now apply the results to the high latitude Ionospheric E region. We 

choose v /$) - 10 v /0 * 10 ^ and V. ' 20 c . For V. - 20 c_ the 
11 ee ’ d s a s 

corresponding shear length Is about 1000 km (N0 + Ions, n - IQ 5 cm -3 ). We 

find that the growth rate Is 0.9 s - *, for modes with wavelengths of 50 

meters (k p - 0.3, and p - 2.5 meters), which Is comparable with the local 

y s s 

growth rate. From the plots of the wave functions we see that the mode Is 

localized In the north-south direction In a region of Ax ** 200 p g * 500 m. 

Moreover, we find that k » k (k - 1/Ax) which Indicates the strong two 

y x x 

dimensional structure of the mode in the plane perpendicular to the magnetic 
field during the linear phase of the instability. 

We conclude that the magnetic shear corresponding to normal auroral 

2 

conditions with field aligned currents of the order of pA/m is small 

(L g ' 1000 km) and thus does not apprecibly alter the current driven 

collisional ion cyclotron mode discussed by Chaturvedi (1976). Even strong 

2 

currents under some disturbed conditions, *» 135 pA/m , reported by Burke et 
al. (1983) do not produce significant shear (L g - 300 km) in the magnetic 
field to have any effect on the CIC. Currents of the order mA/m 2 or more 
possibly produce significant effects. 

Finally, we discuss the limitations of the present theory. (1) The 

mode structure equation (Eq. 18) and the subsequent dispersion relation (Eq. 

23) were derived under the assumption that k z /k y < 1 based on the premise 

that the ion parallel motion term in the ion continuity equation] 

may not be important. Thus the results presented do not contain ion 

parallel motion effects. However, we extended the theory to include the ion 

parallel motion and found no significant changes in the way the magnetic 

shear affects the mode. The local growth rate (in the limit p /L ** 0 ) 

s s 

with the ion parallel motion is smaller by about 10%; the growth rate with 

magnetic shear is consistently smaller but has similar behaviour as shown 

earlier. (2) For larger shears, p /L '0.1, the parabolic expansion of the 

s s 

potential function (Eq. 23), used to derive the analytical dispersion 
relation, may not be adequate. This is due to the fact that the wave 
function spreads out and samples non-parabolic part of the potential. 
Numerical solutions of (Eq. 18) do confirm the spreading of the wave 
packet. Furthermore, the wave packet localizes in the region where the 
effective (k z /ky) '1 because magnetic shear introduces an effective k z * (3) 

9 






We have done an MHD analysis and arrived at the mode structure equation with 
a potential term, Q(x) [Eq. (18)], whose derivative vanishes at x Q , whereas, 
in the collisionless case the derivative of Q(x) at x Q is finite, thus 
drastically affecting the collisionless kinetic ion cyclotron instability 
(Ganguli and Bakshi, 1982). To further examine the above three aspects we 
will present the kinetic analysis of the collisonal ion cyclotron waves in 
the presence of magnetic shear in a future paper. For higher drift 
velocities (comparable to the E * B drift velocities) where the self- 
consistent magnetic shear plays a significant role in determining the mode 
structure, a complete treatment is needed which includes an electric field 
along the y-direction and the related velocity shear. 

In conclusion, we have examined the influence of magnetic shear on the 
collisional current driven ion cyclotron instability. Self consistent 
magnetic shear corresponding to moderate drift velocities (V^ 0.5 - 5 

km/s) near the threshold velocity for the instability does not have 
significant effects on the instability. We find that the mode is almost 
perpendicularly propagating and is localized in a region extending upto a 
few hundred kilometers in the direction of the magnetic field under auroral 

conditions. However, in domains where the plama density is high such that 
2 2 

the parameter a, pe / lil ce is large, we find that the growth rate of a 
particular mode maximizes at an optimum drift velocity much larger than the 
sound speed. Finally, we find that strong shears (p s /L g > 0.05) signifi¬ 
cantly reduce the growth rate and stabilize the collisional current driven 
ion cyclotron instability. This strong shear corresponds to large parallel 
drift velocities (V^ » c g ), as seen in fig. 4. 
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Figure 1. Plot of normalized growth rate (y/0^) versus the normalized • 

wave number (kyP 8 )* The parameters used are v e A* e * 0.01, 
v i /n i ■ 0.01, and v ( j/ c s '* 50. The solid line represents the 
shear free case, and the dashed line represents the case where 

P a /L a - 0.01. • 
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Figure 2. Normalized growth rate (y/Q^) versus P g /L g treating P g /^ 8 as an 

external parameter for V^/Cg ■ 50, k ^p g ■ 0.3, and 0.5, and for 

v /tJ - 0.01, and v./fl, - 0.01. 
e e 11 
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Figure 4. Growth rate versus the drift velocity curves with and without 

self-consistent magnetic shear for the following parameters: 

k p « 0.3, V /0 - 0.01, and v./n, - 0.01. Curve A represents 

y s e e gi l 

the shear-free case and a m 10 . Curves B and C are for a ■ 

10 3 and 10‘\ respectively, and these correspond to /<■)? ■ 

. 1 pe ce 

10 3 and 10 . 
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OICY ATTN J-3 WWMCCS EVALUATION OFFICE 

DIRECTOR 

JOINT STRAT TGT PLANNING STAFF 
OFFUTT AFB 
OMAHA, NB 68113 
OICY ATTN JLTW-2 
OICY ATTN JPST G. GOETZ 

CHIEF 

LIVERMORE DIVISION FLD COMMAND DNA 
DEPARTMENT OF DEFENSE 
LAWRENCE LIVERMORE LABORATORY 
P.O. BOX 808 
LIVERMORE, CA 94550 
OICY ATTN FCPRL 

COMMANDANT 

NATO SCHOOL (SHAPE) 

APO NEW YORK 09172 

OICY ATTN U.S. DOCUMENTS OFFICER 

UNDER SECY OF DEF FOR RSCH 6 ENGRG 
DEPARTMENT OF DEFENSE 
WASHINGTON, D.C. 20301 

OICY ATTN STRATEGIC 6 SPACE SYSTEMS (OS) 

WWMCCS SYSTEM ENGINEERING ORG 
WASHINGTON, D.C. 20305 
OICY ATTN R. CRAWFORD 

COMMANDER/DIRECTOR 
ATMOSPHERIC SCIENCES LABORATORY 
U.S. ARMY ELECTRONICS COMMAND 
WHITE SANDS MISSILE RANGE, NM 88002 
OICY ATTN DELAS-EO F. NILES 


COMMANDER 
FIELD COMMAND 
DEFENSE NUCLEAR AGENCY 
KIRTLAND, AFB, NM 87115 
OICY ATTN FCPR 

DEFENSE NUCLEAR AGENCY 
SAO/DNA 
BUILDING 20676 
KIRTLAND AFB, NM 87115 
OICY D.C. THORNBURG 
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DIRECTOR 

HMD ADVANCED TECH CTR 
HUNTSVILLE OFFICE 
P.O. BOX 1500 
HUNTSVILLE, AL 35807 

01CY ATTN ATC-T MELVIN T. CAPPS 

OICY ATTN ATC-0 W. DAVIES 

OICY ATTN ATC-R DON RUSS 

PROGRAM MANAGER 
BMD PROGRAM OFFICE 
5001 EISENHOWER AVENUE 
ALEXANDRIA, VA 22333 

01CY ATTN DACS-BMT J. SHEA 

CHIEF C-E- SERVICES DIVISION 
U.S. ARMY COMMUNICATIONS CMD 
PENTAGON RM 1B269 
WASHINGTON, D.C. 20310 

01CY ATTN C- E-SERVICES DIVISION 

COMMANDER 

FRADCOM TECHNICAL SUPPORT ACTIVITY 
DEPARTMENT OF THE ARMY 
FORT MONMOUTH, N.J. 07703 

01CY ATTN DRSEL-NL-RD H. BENNET 
01CY ATTN DRSEL-PL-ENV H. BOM ICE 
OICY ATTN J.E. QUIGLEY 

COMMANDER 

U.S. ARMY COMM-ELEC ENGRG INSTAL AGY 
FT. HUACHUCA, AZ 85613 

OICY ATTN CCC-EMEO GEORGE LANE 

COMMANDER 

U.S. ARMY FOREIGN SCIENCE & TECH CTR 
220 7TH STREET, NE 
CHARLOTTESVILLE, VA 22901 
OICY ATTN DRXST-SD 

COMMANDER 

U.S. ARMY MATERIAL DEV & READINESS CMD 
5001 EISENHOWER AVENUE 
ALEXANDRIA, VA 22333 

OICY ATTN DRCLDC J.A. BENDER 

COMMANDER 

U.S. ARMY NUCLEAR AND CHEMICAL AGENCY 
7500 BACKLICK ROAD 
BLDG 2073 

SPRINGFIELD, VA 22150 
OICY ATTN LIBRARY 


DIRECTOR 

U.S. ARMY BALLISTIC RESEARCH LABORATORY 
ABERDEEN PROVING GROUND, MD 21005 
OICY ATTN TECH LIBRARY EDWARD BAICY 

COMMANDER 

U.S. ARMY SATCOM AGENCY 
FT. MONMOUTH, NJ 07703 

OICY ATTN DOCUMENT CONTROL 

COMMANDER 

U.S. ARMY MISSILE INTELLIGENCE AGENCY 
REDSTONE ARSENAL, AL 35809 
OICY ATTN JIM GAMBLE 

DIRECTOR 

U.S. ARMY TRADOC SYSTEMS ANALYSIS ACTIVITY 
WHITE SANDS MISSILE RANGE, NM 88002 
OICY ATTN ATAA-SA 
OICY ATTN TCC/F. PAYAN JR. 

OICY ATTN ATTA-TAC LTC J. HESSE 

COMMANDER 

NAVAL ELECTRONIC SYSTEMS COMMAND 
WASHINGTON, D.C. 20360 

OICY ATTN NAVALEX 034 T. HUGHES 
OICY ATTN PME 117 
OICY ATTN PME 117-T 
OICY ATTN CODE 5011 

COMMANDING OFFICER 
NAVAL INTELLIGENCE SUPPORT CTR 
4301 SUITLAND ROAD, BLDG. 5 
WASHINGTON, D.C. 20390 

OICY ATTN MR. DUBBIN STIC 12 

OICY ATTN NISC-50 

OICY ATTN CODE 5404 J. GALET 

COMMANDER 

NAVAL OCCEAN SYSTEMS CENTER 
SAN DIECO, CA 92152 
OICY ATTN J. FERGUSON 


18 



SAUL RESEARCH LABORATORY 

WASHINGTON, D.C. 20375 

OICY ATTN CODE 4700 S. L. Ossakow 

26 CYS IF UNCLASS- 1 CY IF CLASS) 
OICY ATTN CODE 4701 I Vltkovitsky 
OICY ATTN CODE 4780 J. Huba (100 
CYS IF UNCLASS, I CY IF CLASS) 
OICY ATTN CODE 7500 

OICY ATTN CODE 7550 

OICY ATTN CODE 7580 

OICY ATTN CODE 7551 

OICY ATTN CODE 7555 

OICY ATTN CODE 4730 E. MCLEAN 

OICY ATTN CODE 4108 

OICY ATTN CODE 4730 B. RIPIN 

20CY ATTN CODE 2628 

COMMANDER 

NAVAL SEA SYSTEMS COMMAND 

WASHINGTON, D.C. 20362 
OICY ATTN CAPT R. PITKIN 

COMMANDER 

NAVAL SPACE SURVEILLANCE SYSTEM 

DAHLGREN, VA 22448 

OICY ATTN CAPT J.H. BURTON 


COMMANDER 

AEROSPACE DEFENSE COMMAND/DC 
DEPARTMENT OF THE AIR FORCE 
ENT AFB, CO 80912 

OICY ATTN DC MR. LONG 

COMMANDER 

AEROSPACE DEFENSE COMMAND/XPD 
DEPARTMENT OF THE AIR FORCE 
ENT AFB, CO 80912 

OICY ATTN XPDQQ 

OICY ATTN XP 

AIR FORCE GEOPHYSICS LABORATORY 
HANSCOM AFB, MA 01731 

OICY ATTN OPR HAROLD GARDNER 
OICY ATTN LKB KENNETH S.W. CHAMPION 
OICY ATTN OPR ALVA T. STAIR 

OICY ATTN PHD JURGEN BUCHAU 

OICY ATTN PHD JOHN P. MULLEN 

AF WEAPONS LABORATORY 
KIRTLAND AFT, NM 87117 
OICY ATTN SUL 

OICY ATTN CA ARTHUR H. GUENTHER 
OICY ATTN NTYCE 1LT. G. KRAJEI 



OFFICER-IN-CHARGE 
NAVAL SURFACE WEAPONS CENTER 
WHITE OAK, SILVER SPRING. MD 20910 
OICY ATTN CODE F31 

DIRECTOR 

STRATEGIC SYSTEMS PROJECT OFFICE 
DEPARTMENT OF THE NAVY 
WASHINGTON, D.C. 20376 
OICY ATTN NSP-2141 
OICY ATTN NSSP-2722 FRED WIMBERLY 

COMMANDER 

NAVAL SURFACE WEAPONS CENTER 
DAHLGREN LABORATORY 
DAHLGREN, VA 22448 

OICY ATTN CODE DF-14 R. BUTLER 

OFFICER OF NAVAL RESEARCH 
ARLINGTON, VA 22217 
OICY ATTN CODE 465 
OICY ATTN CODE 461 
OICY ATTN CODE 402 
OICY ATTN CODE 420 
OICY ATTN CODE 421 


AFTAC 

PATRICK AFB, FL 32925 
OICY ATTN TF/MAJ WILEY 
OICY ATTN TN 

AIR FORCE AVIONICS LABORATORY 
WRIGHT-PATTERSON AFB, OH 45433 
OICY ATTN AAD WADE HUNT 
OICY ATTN AAD ALLEN JOHNSON 

DEPUTY CHIEF OF STAFF 
RESEARCH, DEVELOPMENT, & ACQ 
DEPARTMENT OF THE AIR FORCE 
WASHINGTON, D.C. 20330 
OICY ATTN AFRDQ 

HEADQUARTERS 

ELECTRONIC SYSTEMS DIVISION 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01731 
OICY ATTN J. DEAS 

HEADQUARTERS 

ELECTRONIC SYSTEMS DIVISION/YSEA 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01732 
OICY ATTN YSEA 
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HEADQUARTERS 

ELECTRONIC SYSTEMS DIVISION/DC 
DEPARTMENT OF THE AIR FORCE 
HANSCOM AFB, MA 01731 

01CY ATTN DCKC MAJ J.C. CLARK 

COMMANDER 

FOREIGN TECHNOLOGY DIVISION, AFSC 
WRIGHT-PATTERSON AFB, OH 45433 
01CY ATTN NICD LIBRARY 

01CY ATTN ETDP B. BALLARD 

COMMANDER 

ROME AIR DEVELOPMENT CENTER, AFSC 
GRIFFISS AFB, NY 13441 

01CY ATTN DOC LIBRARY/TSLD 
01CY ATTN OCSE V. COYNE 

SAMSO/SZ 

POST OFFICE BOX 92960 
WORLDUAY POSTAL CENTER 
LOS ANGELES, CA 90009 
(SPACE DEFENSE SYSTEMS) 

01CY ATTN SZJ 

STRATEGIC AIR COMMAND/XPFS 
OFFUTT AFB, NB 68113 

01CY ATTN ADWATE MAJ BRUCE BAUER 
01CY ATTN NRT 

01CY ATTN DOK CHIEF SCIENTIST 

SAMSO/SK 
P.O. BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 

01CY ATTN SKA (SPACE COMM SYSTEMS) 
M. CLAVIN 

SAMSO/MN 

NORTON AFB, CA 92409 
(MINUTEMAN) 

01CY ATTN MNNL 

COMMANDER 

ROME AIR DEVELOPMENT CENTER, AFSC 
HANSCOM AFB, MA 01731 

01CY ATTN EEP A. LORENTZEN 

DEPARTMENT OF ENERGY 
LIBRARY ROOM G-042 
WASHINGTON, D.C. 20545 

01CY ATTN DOC CON FOR A. LABOWITZ 


DEPARTMENT OF ENERGY 
ALBUQUERQUE OPERATIONS OFFICE 
P.O. BOX 5400 
ALBUQUERQUE, NM 87115 

01CY ATTN DOC CON FOR D. SHERWOOD 

EG&G, INC. 

LOS ALAMOS DIVISION 

P.O. BOX 809 

LOS ALAMOS, NM 85544 

01CY ATTN DOC CON FOR J. BREEDLOVE 

UNIVERSITY OF CALIFORNIA 
LAWRENCE LIVERMORE LABORATORY 
P.O. BOX 808 
LIVERMORE, CA 94550 

01CY ATTN DOC CON FOR TECH INFO DEPT 

01CY ATTN DOC CON FOR L-389 R. OTT 

01CY ATTN DOC CON FOR L-31 R. HAGER 

01CY ATTN DOC CON FOR L-46 F. SEWARD 

LOS ALAMOS NATIONAL LABORATORY 

P.O. BOX 1663 

LOS ALAMOS, NM 87545 

01CY ATTN DOC CON FOR J. WOLCOTT 

01CY ATTN DOC CON FOR R.F. TASCHEK 

01CY ATTN DOC CON FOR E. JONES 

01CY ATTN DOC CON FOR J. MALIK 

01CY ATTN DOC CON FOR R. JEFFRIES 

01CY ATTN DOC CON FOR J. ZINN 

01CY ATTN DOC CON FOR P. KEATON 

01CY ATTN DOC CON FOR D. WESTERVELT 
01CY ATTN D. SAPPENFIELD 

SAND1A LABORATORIES 
P.O. BOX 5800 
ALBUQUERQUE, NM 87115 

01CY ATTN DOC CON FOR W. BROWN 

OICY ATTN DOC CON FOR A. THORNBROUGH 

OICY ATTN DOC CON FOR T. WRIGHT 

OICY ATTN DOC CON FOR D. DAHLGREN 

OICY ATTN DOC CON FOR 3141 

OICY ATTN DOC CON FOR SPACE PROJECT DIV 

SANDIA LABORATORIES 
LIVERMORE LABORATORY 
P.O. BOX 969 
LIVERMORE, CA 94550 

OICY ATTN DOC CON FOR B. MURPHEY 

OICY ATTN DOC CON FOR T. COOK 


OFFICE OF MILITARY APPLICATION 
DEPARTMENT OF ENERGY 
WASHINGTON, D.C. 20545 

OICY ATTN DOC CON DR. YO SONG 
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OTHER GOVERNMENT 

DEPARTMENT OF COMMERCE 
NATIONAL BUREAU OF STANDARDS 
WASHINGTON, D.C. 20234 

OICY (ALL CORRES: ATTN SEC OFFICER FOR) 

INSTITUTE FOR TELECOM SCIENCES 
NATIONAL TELECOMMUNICATIONS & INFO ADMIN 
BOULDER, CO 80303 

01CY ATTN A. JEAN (UNCLASS ONLY) 

01CY ATTN W. UTLAUT 
01CY ATTN D. CROMBIE 
01CY ATTN L. BERRY 

NATIONAL OCEANIC & ATMOSPHERIC ADMIN 
ENVIRONMENTAL RESEARCH LABORATORIES 
DEPARTMENT OF COMMERCE 
BOULDER, CO 80302 
01CY ATTN R. GRUBB 
01CY ATTN AERONOMY LAB G. REID 


DEPARTMENT OF DEFENSE CONTRACTORS 

AEROSPACE CORPORATION 
P.O. BOX 92957 
LOS ANGELES, CA 90009 
01CY ATTN I. GARFUNKEL 
01CY ATTN T. SALMI 
01CY ATTN V. JOSEPHSON 
01CY ATTN S. BOWER 
01CY ATTN D. OLSEN 

ANALYTICAL SYSTEMS ENGINEERING CORP 
5 OLD CONCORD ROAD 
BURLINGTON, MA 01803 

01CY ATTN RADIO SCIENCES 

AUSTIN RESEARCH ASSOC., INC. 

1901 RUTLAND DRIVE 
AUSTIN, TX 78758 
01CY ATTN L. SLOAN 
OICY ATTN R. THOMPSON 

BERKELEY RESEARCH ASSOCIATES, INC. 
P.O. BOX 983 
BERKELEY, CA 94701 
OICY ATTN J. WORKMAN 
OICY ATTN C. PRETTIE 
OICY ATTN S. BRECHT 


BOEING COMPANY, THE 
P.O. BOX 3707 
SEATTLE, WA 98124 
OICY ATTN G. KEISTER 
OICY ATTN D. MURRAY 
OICY ATTN G. HALL 
OICY ATTN J. KENNEY 

CHARLES STARK DRAPER LABORATORY, INC. 
555 TECHNOLOGY SQUARE 
CAMBRIDGE, MA 02139 
OICY ATTN D.B. COX 
OICY ATTN J.P. GILMORE 

COMSAT LABORATORIES 
LINTHICUM ROAD 
CLARKSBURG, MD 20734 
OICY ATTN G. HYDE 

CORNELL UNIVERSITY 

DEPARTMENT OF ELECTRICAL ENGINEERING 
ITHACA, NY 14850 

OICY ATTN D.T. FARLEY, JR. 

ELECTROSPACE SYSTEMS, INC. 

BOX 1359 

RICHARDSON, TX 75080 
OICY ATTN H. LOGSTON 
OICY ATTN SECURITY (PAUL PHILLIPS) 

EOS TECHNOLOGIES, INC. 

606 Ullshlre Blvd. 

Santa Monica, Calif 90401 
OICY ATTN C.S. GABBARD 

ESL, INC. 

495 JAVA DRIVE 
SUNNYVALE, CA 94086 
OICY ATTN J. ROBERTS 
OICY ATTN JAMES MARSHALL 

GENERAL ELECTRIC COMPANY 
SPACE DIVISION 
VALLEY FORGE SPACE CENTER 
GODDARD BLVD KING OF PRUSSIA 
P.O. BOX 8555 
PHILADELPHIA, PA 19101 

OICY ATTN M.H. BORTNER SPACE SCI LAB 

GENERAL ELECTRIC COMPANY 
P.O. BOX 1122 
SYRACUSE, NY 13201 
OICY ATTN F. REIBERT 
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general electric TECH SERVICES CO. INC. 
RMES 

COl'RT STREET 
SYRACUSE, NY 13201 
01CY ATTN G. MILLMAN 

GEOPHYSICAL INSTITUTE 
UNIVERSITY OF ALASKA 
FAIRBANKS, AK 99701 

(ALL CLASS ATTN: SECURITY OFFICER) 
OICY ATTN T.N. DAVIS (UNCLASS ONLY) 
01CY ATTN TECHNICAL LIBRARY 
OICY ATTN NEAL BROWN (UNCLASS ONLY) 

CTE SYLVANIA, INC. 

ELECTRONICS SYSTEMS CRP-EASTERN DIV 
77 A STREET 
NEEDHAM, MA 02194 

OICY ATTN DICK STEINHOF 

HSS, INC. 

2 ALFRED CIRCLE 
BEDFORD, MA 01730 

OICY ATTN DONALD HANSEN 

ILLINOIS, UNIVERSITY OF 
107 COBLE HALL 
150 DAVENPORT HOUSE 
CHAMPAIGN, IL 61820 

(ALL CORRES ATTN DAN MCCLELLAND) 

OICY ATTN K. YEH 

INSTITUTE FOR DEFENSE ANALYSES 
1801 NO. BEAUREGARD STREET 
ALEXANDRIA, VA 22311 
OICY ATTN J.M. AEIN 
OICY ATTN ERNEST BAUER 
OICY ATTN HANS WOLFARD 
OICY ATTN JOEL BENGSTON 

INTL TEL 6 TELEGRAPH CORPORATION 
500 WASHINGTON AVENUE 
NUTLEY, NJ 07110 
OICY ATTN TECHNICAL LIBRARY 

JAYCOR 

11011 TORREYANA ROAD 
P.O. BOX 85154 
SAN DIEGO, CA 92138 
OICY ATTN J.L. SPERLING 


JOHNS HOPKINS UNIVERSITY 
APPLIED PHYSICS LABORATORY 
JOHNS HOPKINS ROAD 
LAUREL, MD 20810 

OICY ATTN DOCUMENT LIBRARIAN 
OICY ATTN THOMAS POTEMRA 
OICY ATTN JOHN DASSOULAS 

KAMAN SCIENCES CORP 
P.O. BOX 7463 

COLORADO SPRINGS, CO 80933 
OICY ATTN- T. MEAGHER 

KAMAN TEMPO -CENTER FOR ADVANCED STUDIES 
816 STATE STREET (P.O DRAWER QQ) 

SANTA BARBARA, CA 93102 
OICY ATTN DASIAC 
OICY ATTN WARREN S. KNAPP 
OICY ATTN WILLIAM MCNAMARA 
OICY ATTN B. GAMBILL 

LINKABIT CORP 
10453 ROSELLE 
SAN DIEGO, CA 92121 
OICY ATTN IRWIN JACOBS 

LOCKHEED MISSILES & SPACE CO., INC 
P.O. BOX 504 
SUNNYVALE, CA 94088 
OICY ATTN DEPT 60-12 
OICY ATTN D.R. CHURCHILL 

LOCKHEED MISSILES & SPACE CO., INC. 

3251 HANOVER STREET 
PALO ALTO, CA 94304 

OICY ATTN MARTIN WALT DEPT 52-12 
OICY ATTN W.L. IMHOF DEPT 52-12 
OICY ATTN RICHARD G. JOHNSON DEPT 52-12 
OICY ATTN J.B. CLADIS DEPT 52-12 

MARTIN MARIETTA CORP 
ORLANDO DIVISION 
P.O. BOX 5837 
ORLANDO, FL 32805 
OICY ATTN R. HEFFNER 


M.I.T. LINCOLN LABORATORY 
P.O. BOX 73 
LEXINGTON, MA 02173 

OICY ATTN DAVID M. TOWLE 
OICY ATTN L. LOUGHLIN 
OICY ATTN D. CLARK 
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MCDONNEL DOUGLAS CORPORATION 
5301 BOLSA AVENUE 
HUNTINGTON BEACH, CA 92647 
OICY ATTN N. HARRIS 
OICY ATTN J. MOULE 
OICY ATTN GEORGE MROZ 
OICY ATTN W. OLSON 
OICY ATTN R.W. HALPRIN 
OICY ATTN TECHNICAL LIBRARY SERVICES 

MISSION RESEARCH CORPORATION 
735 STATE STREET 
SANTA BARBARA, CA 93101 
OICY ATTN P. FISCHER 
OICY ATTN W.F. CREVIER 
OICY ATTN STEVEN L. GUTSCHE 
OICY ATTN R. BOGUSCH 
OICY ATTN R. HENDRICK 
OICY ATTN RALPH KILB 
OICY ATTN DAVE SOWLE 
OICY ATTN F. FAJEN 
OICY ATTN M. SCHEIBE 
OICY ATTN CONRAD L. LONGMIRE 
OICY ATTN B. WHITE 

MISSION RESEARCH CORP. 

1720 RANDOLPH ROAD, S.E. 

ALBUQUERQUE, NEW MEXICO 87106 
OICY R. STELLINGWERF 
OICY M. ALME 
OICY L. WRIGHT 

MITRE CORPORATION, THE 
P.O. BOX 208 
BEDFORD, MA 01730 

OICY ATTN JOHN MORGANSTERN 
OICY ATTN G. HARDING 
OICY ATTN C.E. CALLAHAN 

MITRE CORP 

WESTGATE RESEARCH PARK 
1820 DOLLY MADISON BLVD 
MCLEAN, VA 22101 
OICY ATTN W. HALL 
OICY ATTN W. FOSTER 

PACIFIC-SIERRA RESEARCH CORP 
12340 SANTA MONICA BLVD. 

LOS ANGELES, CA 90025 

OICY ATTN E.C. FIELD, JR. 


PENNSYLVANIA STATE UNIVERSITY 
IONOSPHERE RESEARCH LAB 
318 ELECTRICAL ENGINEERING EAST 
UNIVERSITY PARK, PA 16802 
(NO CLASS TO THIS ADDRESS) 

OICY ATTN IONOSPHERIC RESEARCH LAB 

PHOTOMETRICS, INC. 

4 ARROW DRIVE 
WOBURN, MA 01801 

OICY ATTN IRVING L. KOFSKY 

PHYSICAL DYNAMICS, INC. 

P.O. BOX 3027 
BELLEVUE, WA 98009 

OICY ATTN E.J. FREMOUW 

PHYSICAL DYNAMICS, INC. 

P.O. BOX 10367 
OAKLAND, CA 94610 
ATTN A. THOMSON 

R 6 D ASSOCIATES 

P.O. BOX 9695 

MARINA DEL REY, CA 90291 


OICY 

ATTN 

FORREST GILMORE 

OICY 

ATTN 

WILLIAM B. WRIGHT, 

OICY 

ATTN 

ROBERT F. LELEVIER 

OICY 

ATTN 

WILLIAM J. KARZAS 

OICY 

ATTN 

H. 

ORY 

OICY 

ATTN 

C. 

MACDONALD 

OICY 

ATTN 

R. 

TURCO 

OICY 

ATTN 

L. 

DeRAND 

OICY 

ATTN 

W. 

TSAI 


RAND CORPORATION, THE 
1700 MAIN STREET 
SANTA MONICA, CA 90406 
OICY ATTN CULLEN CRAIN 
OICY ATTN ED BEDROZIAN 

RAYTHEON CO. 

528 BOSTON POST ROAD 
SUDBURY, MA 01776 

OICY ATTN BARBARA ADAMS 

RIVERSIDE RESEARCH INSTITUTE 
330 WEST 42nd STREET 
NEW YORK, NY 10036 

OICY ATTN VINCE TRAPANI 
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SCIENCE APPLICATIONS, INC. 
1150 PROSPECT PLAZA 
LA JOLLA, CA 92037 

OICY ATTN LEWIS M. LINSON 
OICY ATTN DANIEL A. HAMLIN 
01CY ATTN E. FRIEMAN 
01CY ATTN E.A. STRAKER 
OICY ATTN CURTIS A. SMITH 
OICY ATTN JACK MCDOUGALL 

SCIENCE APPLICATIONS, INC 
1710 GOODRIDGE DR. 

MCLEAN, VA 22102 
ATTN: J. COCKAYNE 

SRI INTERNATIONAL 
333 RAVENSWOOD AVENUE 
MENLO PARK, CA 94025 

OICY ATTN DONALD NEILSON 
OICY ATTN ALAN BURNS 
OICY ATTN G. SMITH 
OICY ATTN R. TSUNODA 
OICY ATTN DAVID A. JOHNSON 
OICY ATTN WALTER G. CHESNUT 
OICY ATTN CHARLES L. RINO 
OICY ATTN WALTER JAYE 
OICY ATTN J. VICKREY 
OICY ATTN RAY L. LEADABRAND 
OICY ATTN G. CARPENTER 
OICY ATTN G. PRICE 
OICY ATTN R. LIVINGSTON 
OICY ATTN V. GONZALES 
OICY ATTN D. MCDANIEL 

TECHNOLOGY INTERNATIONAL CORP 
75 WIGGINS AVENUE 
BEDFORD, MA 01730 
OICY ATTN W.P. BOQUIST 

TOYON RESEARCH CO. 

P.O. Box 6890 
SANTA BARBARA, CA 93111 
OICY ATTN JOHN ISE, JR. 

OICY ATTN JOEL GARBARINO 


TRW DEFENSE & SPACE SYS GROUP 
ONE SPACE PARK 
REDONDO BEACH, CA 90278 
OICY ATTN R. K. PLEBUCH 
OICY ATTN S. ALTSCHULER 
OICY ATTN D. DEE 
OICY ATTN D/ STOCKWELL 
SNTF/1575 

VI SI DYNE 

SOUTH BEDFORD STREET 
BURLINGTON, MASS 01803 
OICY ATTN W. REIDY 
OICY ATTN J. CARPENTER 
OICY ATTN C. HUMPHREY 
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IONOSPHERIC MODELING DISTRIBUTION LIST 
(UNCLASSIFIED ONLY) 

PLEASE DISTRIBUTE ONE COPY TO EACH OF THE FOLLOWING PEOPLE (UriLESS OTHERWISE 
NOTED) 


NAVAL RESEARCH LABORATORY 
WASHINGTON, D.C. 20375 

Dr. P. MANGE - CODE 4101 
Dr. P. GOODMAN - CODE 4180 

A.F. GEOPHYSICS LABORATORY 
L.G. HANSCOM FIELD 
BEDFORD, MA 01730 
DR. T. ELKINS 
DR. W. SWIDER 
MRS. R. SAGALYN 
DR. J.M. FORBES 
DR. T.J. KENESHEA 
DR. W. BURKE 
DR. H. CARLSON 
DR. J. JASPERSE 

BOSTON UNIVERSITY 
DEPARTMENT OF ASTRONOMY 
BOSTON, MA 02215 
DR. J. AARONS 

CORNELL UNIVERSITY 
ITHACA, NY 14850 
DR. W.E. SWARTZ 
DR. D. FARLEY 
DR. M. KELLEY 

HARVARD UNIVERSITY 
HARVARD SQUARE 
CAMBRIDGE, MA 02138 
DR. M.B. HcELROY 
DR. R. LINDZEN 


NASA 

GODDARD SPACE FLICHT CENTER 
GREENBELT, MD 20771 
DR. K. MAEDA 
DR. S. CURTIS 
DR. M. DUBIN 

DR. N. MAYNARD - CODE 696 
COMMANDER 

NAVAL AIR SYSTEMS COMMAND 
DEPARTMENT OF THE NAVY 
WASHINGTON, D.C. 20360 
DR. T. CZUBA 

COMMANDER 

NAVAL OCEAN SYSTEMS CENTER 
SAN DIEGO, CA 92152 

MR. R. ROSE - CODE 5321 

NOAA 

DIRECTOR OF SPACE AND 
ENVIRONMENTAL LABORATORY 
BOULDER, CO 80302 
DR. A. GLENN JEAN 
DR. G.W. ADAMS 
DR. D.N. ANDERSON 
DR. K. DAVIES 
DR. R.F. DONNELLY 

OFFICE OF NAVAL RESEARCH 
800 NORTH QUINCY STREET 
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